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Flip-flop phenomenon: observations and theory 
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Abstract. In many active stars the spots concentrate on two permanent active longitudes which are 180° apart. In some of 
these stars the dominant part of the spot activity changes the longitude every few years. This so-called flip-flop phenomenon 
has up to now been reported in 1 1 stars, both single and binary alike, and including also the Sun. To explain this phenomenon, 
a non-axisymmetric dynamo mode, giving rise to two permanent active longitudes at opposite stellar hemispheres, is needed 
together with an oscillating axisymmetric magnetic field. Here we discuss the observed characteristics of the flip-flop phe- 
nomenon and present a dynamo solution to explain them. 
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1. Introduction 

The global structure and behaviour of the stellar magnetic 
fields are determined by different dynamo modes that have 
different symmetries and stabilities (see e.g. Brandenburg et 
al.[l989). In slowly rotating stars, like the Sun, axisymmet- 
ric modes are excited. These modes do not show any struc- 
ture in the longitudinal distribution of the spots and they os- 
cillate in time. In more rapidly rotating stars the higher or- 
der non-axisymmetric modes get excited (see e.g. Moss et 
al. 119951 Tuominen, Berdyugina & Korpi 2002). The mag- 
netic configuration in the non-axisymmetric modes consists 
of two starspots that are 180° apart, explaining the perma- 
nent active longitudes seen in many rapidly rotating stars. 
These non-axisymmetric modes do not oscillate. For explain- 
ing the flip-flop phenomenon, where we see both active lon- 
gitudes and oscillations, axisymmetric dynamo modes need 
to co-exist with the non-axisymmetric modes. 

In this paper we describe the observed characteristics of 
the flip-flop phenomenon and present a model that can pro- 
duce them. 

2. Observations of flip-flop phenomenon 

The flip-flop phenomenon, in which the main part of the spot 
activity changes 1 80° on the stellar surface, was first discov- 
ered in the early 1990s on a single, very active, giant, FK 
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Com. Jetsu et al. (119911 119931) noticed from photometric ob- 
servation, that the spot activity on FK Com for 1966-1990 
concentrated on two longitudes, 180° apart. They also noted 
that during the individual observing seasons only one of the 
active longitudes had spots. This behaviour is very well il- 
lustrated in Fig. [I] (from Jetsu et a; 1993), which shows the 
normalized magnitudes of FK Com for 1966-1990. The pho- 
tometric minimum is always either around the phase 0.0 or 
0.5. The data Jetsu et al. (119931) used can be analysed together 
with more recent observations (1991-2003) to estimate the 
frequency at which flip-flop events occur on FK Com. There 
is on average one flip-flop event every 2.6 years, giving full 
cycle length of 5.2 years (Korhonen et al. 2004). 

After the discovery of the flip-flop phenomenon on 
FK Com it has been reported also on other active stars. 
Berdyugina & Tuominen ( 1998) studied the photometric ob- 
servations of four RS CVn binaries and discovered that also 
these stars have permanent active longitudes that are alterna- 
tively active. In the case of II Peg Rodono et al. (2000) later 
confirmed their results. Berdyugina, Pelt & Tuominen ( 2002 ) 
discovered flip-flops on young solar type star, LQ Hya. And 
a recent analysis of 120 years of sunspot data (Berdyugina & 
Usoskin[2003) suggests that the Sun also has permanent ac- 
tive longitudes with associated flip-flops. On the Sun a flip- 
flop event occurs on average every 3.8 years on the northern 
and every 3.65 years on the southern hemisphere. 

When the flip-flop phenomenon was discovered, it was 
not sure whether the phenomenon was caused by spot move- 
ment across the stellar disk or emergence of flux on the new 
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Fig. 1. The active longitude structure on FK Com. Normal- 
ized magnitudes (ordinate) of FK Com for 1966-1990 plot- 
ted against the phase (abscissa). The phases have been deter- 
mined using the ephemeris HJD 2439252.8 95 + 2.4 002466 E. 
This figure has been taken from Jetsu et al. ( 1993 ). 



active longitude. Korhonen et al. (2001) have shown, with 
Doppler images just before and after a flip-flop event on 
FK Com, that flip-flops are caused by changing the relative 
strengths of the spot groups at the two active longitudes with- 
out actual spot movement on the stellar surface. 

All the stars for which the flip-flop phenomenon has been 
reported are listed in Table [l] Spectral type, rotation period 
and the relative differential rotation coefficient are given to- 
gether with the flip-flop period (the length of the full cycle). 
The flip-flop phenomenon has so far been detected in many 
different kinds of stars: both binaries and single stars; young, 
main sequence and evolved alike. Usually, the flip-flop period 
is between 5 and 10 years, median being 7 years. The stars 
themselves usually have rotation periods < 3 days (median 
2.4 days). Anyhow, no clear correlation between the rotation 
period and the flip-flop period can be seen. 

Apart from the stars mentioned in Table [l] there are 
also two other stars for which flip-flops have been reported. 
These stars are single giant HD 199 178 (Hackman 2004| and 
RS CVn binary RT Lac (Lanza et al. 12002ft . In these stars 
only very few events have been observed, so no information 
on the flip-flop cycle length can be obtained. 



Table 1. Stars that show flip-flop phenomenon. In the Table 
the name of the star, spectral type and age, rotation period, 
flip-flop period and the relative differential rotation coeffi- 
cient are given. 



Name 


Type 


Prot 


P ff 




Sun 


single, G2 V 


27 d 


7 V 


0.19 


LQHya 


single, K2 V, ZAMS 


1.6 d 


5.2 2 yr 


0.022 3 


AB Dor 


single, K0 V, ZAMS 


0.5 d 


5.5 4 yr 


0.05 5 


EK Dra 


single, G1.5 V, ZAMS 


2.6 d 


4 4 yr 




FK Com 


single, G7 III 


2.4 d 


5.2 6 yr 


0.018 6 


II Peg 


RS CVn, K2 IV 


6.7 d 


9.3 7 yr 


0.04 8 


sigma Gem 


RS CVn, Kl III 


19.6 d 


14.9 7 yr 


< 0.004 9 


EI Eri 


RS CVn, G5 IV 


1.95 d 


9.0 7 yr 


-0.15--0.20 10 


HR 7275 


RS CVn, Kl III-IV 


2.3 d 


17.5 7 yr 




l)Berdyugina 


& Usoskin 120031 2)Berdyu 


gina et 


al. 120021 


3)K6vari et al. 200^ 



4)Berdyugina & Jarvinen 2005 5)Collier Cameron & Donati 2002 6)Korhonen 
et al. [2004 7)Berdyugina & Tuominen fl998l 8)Weber 120041 9)K6vari et al. I200T1 
10)Washuettl 2004 
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Fig. 2. The solar rotation law. 
3. Modelling flip-flops 

The model consists of a turbulent fluid in a spherical shell of 
inner radius r- in and outer radius r ou t • 
We solve the induction equation 

d(B) 



dt 



curl(a ? o (B) — ryxcurl(B)), 



(1) 



in spherical coordinates (r, (p) for an a 2 ^-dynamo. A solar 
type rotation law (see Fig.|3 in the corotating frame with the 
core 



1 



1 + erf , 

di 

where ft s = Vt eq — acos 2 is used. 



(2) 



a 



Only the symmetric part 

rr = aoCOS#(l. — 2cOS 2 #) 

aee — aocos#(l. — 2sin 2 #) 
a r e = otdr = 2aocos 2 #sin# 



(3) 

of the a-tensor is included. In order to saturate the dynamo 
we choose a local quenching of 

(4) 



1 



B 2 /B 2 



eq 



For ao we choose a value slightly above the critical one for 
the dynamo threshold. 

The inner boundary is a perfect conductor and the outer 
boundary resembles a vacuum condition, by including an 
outer region up to 1 .2 stellar radii into the computational grid 
with 10 times higher diffusivity. At the very outer part the 
pseudo vacuum condition is used. In order to see the influ- 
ence of the thickness of the convection zone we have chosen 
ri n = 0.7 for a thin (results shown in Fig.O and ri n = 0.4 
for a thick (Fig.© convection zone. 



4. Results 

With the parameter Qo we model the strength of the differen- 
tial rotation. For Flo = 1 we have a solar rotation law with an 
oscillating axisymmetric dynamo solution. The case = 
is an a 2 -dynamo which gives a migrating non-axisymmetric 
dynamo because of the anisotropic a (cf . Riidiger, Elstner & 



Ossendrijver 2003 
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Fig. 3. The magnetic pressure from the dominant radial component of the magnetic field on the stellar surface is shown at 
three different time steps for the thin layer model. We subtracted rotation and migration. 




Fig. 4. The upper panel is as in Fig.[3]but now for the thick convection zone model. The lower panel shows Doppler images 
of FK Com for June 1997 and January 1998 (from Korhonen et al. 2001). For the Doppler images the grey scale gives the 
temperature scale of 3600 K - 5700 K. 



For 10% of the solar differential rotation we found similar 
excitation conditions for a drifting non-axisymmetric mode 
and an oscillating axisymmetric mode. Because of the cho- 
sen positive a in the northern hemisphere, we get a pole- 
ward migration of the oscillating mode. The drift of the non- 
axisymmetric mode is opposite to the rotation. 

Using a simple a-quenching, given by Eq.|4l in 3D simu- 
lations we found coexisting solutions for both modes, show- 
ing a magnetic flip-flop phenomenon. We followed the solu- 
tion in our simulation up to 100 diffusion times. There were 
no sign for it being only a temporary phenomenon. The tem- 
poral behaviour of the magnetic energy is shown in Fig.EJ 

For an assumed turbulent diffusivity of about 
10 12 cm 2 s _1 we get a period of about 6 years for the 
thin and 9 years for the thick model. These values are 



mainly determined by the magnetic diffusivity and vary only 
weakly with the value for a. The field strength saturates 
about the equipartition value. In Table 2 we present a choice 
of models in order to illustrate the parameter dependence 
of different solutions for an axisymmetric (first row), a 
non-axisymmetric (second row) and two flip-flop (thick and 
thin model; third and fourth row, respectively) solutions. 
Notice, that the axisymmetric solution appears already for 
^0 = 0.1 but with a higher a than is used for the thick 
convection zone flip-flop model (third row). This is probably 
due to the local quenching of the m = 1 mode. 

For the thick convection zone we found models where 
the magnetic spots appear already at 50° latitude and are mi- 
grating polewards during the cycle. In this case the opposite 
spot does not start to appear exactly 180° away from the old 
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Fig. 5. Energy density in m = and m = 1 modes for the 
thick model normalized to equipartition. The time unit is the 
diffusion time of 30 years. Left: The whole simulation Right: 
Final time with higher resolution. Notice the weak oscillation 
of the energy in the m = 1 mode synchronized with m = 
due to a-quenching. The migration period in Table |2| is di- 
rectly taken from the m = 1 field. 

Table 2. Parameters for a sample of runs with the energy den- 
sity in equipartition units E0 for the mode m = and El for 
m = 1, respectively. P0 denotes the period of the oscillation 
in diffusion times (30 years) and PI the migration period. The 
dynamo-number C a = c^o^star^ -1 is used. 



Tin 


n 


Col 


E0 


El 


P0 


PI 


0.4 


0.1 


20 


1.4 


io- y 


0.27 




0.4 


0.11 


10 


10~ 6 


0.1 




0.23 


0.4 


0.12 


10 


0.06 


0.06 


0.3 


0.2 


0.7 


0.11 


21 


0.1 


0.4 


0.23 


0.46 



spot. The distance can shrink down to 90°. It depends also 
on a. For weakly over critical a we find nearly 180° distance 
between old and new spot. In all models we found a counter- 
rotating migration of the magnetic pattern. The period was 
twice the flip-flop oscillation period for the thin layer and 
nearly equal to the flip-flop period for our thick layer model. 
This migration should be carefully considered in the future 
data analysis. 



4.1. Comparison to FK Com 

Korhonen et al. (2004) found that for FK Com the relative 
surface differential rotation (AO/O) is about 10% of the so- 
lar value. This implies that the 10% scaling of the solar ro- 
tation law, that has been used in the model, is an appropriate 
first approximation for this star. In Fig.|4|the modeling results 
for the thick convection zone model are compared to the ob- 
served surface temperature distributions on FK Com for June 
1997 and January 1998 (Korhonen et al. l2Q0Tl . In the calcula- 
tions the spots appear at high latitudes. This is also consistent 
with the Doppler imaging results for FK Com, which show 
spots mainly at high latitudes, and basically never lower than 
latitude 45°. Also, the spot size looks reasonable and is in 
line with the observations. 



5. Discussion 

For the first time we found stable mixed mode solutions for 
weakly differential rotating stars. We could follow the cycle 
over 100 diffusion times. 

Moss (2004) presented a similar model with isotropic a. 
Because there is no preferred m=l mode in that case, it is 
probably not possible to have stable mixed mode solutions 
for a longer time. Also he did not find similar values for the 
magnetic energy in both modes. In contrary to our solution 
he got an asymmetric distribution of m=l and m = modes 
with respect to the equator. We observed a similar exotic be- 
haviour for highly over-critical a. 

The assumptions for the model are somewhat uncertain. 
First, changing the diffusivity would change the period of the 
oscillating mode and therefore also the flip-flop period. The 
simple scaling of the solar rotation law may not be the best 
approximation for the stars showing flip-flops, also the simple 
a-quenching non-linearity may not be adequate. Anyhow, the 
main properties remain for different quenching forms. Also, 
the simple diagonal form of the a tensor is not justified. Nev- 
ertheless, the results from our model are encouraging. 

The flip-flop phenomenon appears in a limited range of 
the strength of differential rotation. The thickness of the con- 
vection zone is not very important. To what extent a large 
scale meridional flow changes this behaviour has to be inves- 
tigated. 
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